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TRADEMARK OFFICE 
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TITLE: METHOD AND APPARATUS FOR SYNCHRONIZATION OF 
CLOCKS 

INVENTOR: LIN JIN. Beijing, CHINA 

RELATED APPLICATION: 

This application is based on provisional application serial no. 60/271,950 filed on 
03/01/2001 and claims priority of that date. The entire content of the provisional 
application is incorporated herein by reference. 

FIELD OF THE INVENTION: 

This application relates to the synchronization of clocks using, in one application, 
signals from orbiting satellites such as in the Global Positioning System (GPS) and 
GLONASS. 

BACKGROUND OF THE INVENTION: 

In industry there is a business segment which provides electronic devices which 
yield precise time. Some of these devices derive high precision time through a transfer 
process whereby they transfer time from a satellite system. One such satellite system is 
the Global Positioning System (GPS). Various manufacturers use GPS receivers in 
which precise time is derived from one of the atomic clocks in the satellites (GPS Time 
Receiver). These are known by various names such as GPS Time Receivers, Time & 
Frequency Standards, and other names. Sometimes they also provide precise frequency 
measurement. Such Time Receivers are sold for a variety of applications where precise 
time is important. In certain applications there is a need to synchronize these Time 



Docket 02-101 as filed 27FEB2002 



1 



Docket No. 02/101 A 



Receivers with each other. One such application area is the synchronization of 
telecommunications signals. The synchronization process is handled in a number of 
conventional ways. 

When satellite clocks are required to be instantaneously synchronized with each 
other at the nanosecond level difficulties arise due to the misapplication or 
misinterpretation of the relativistic effects and gravity effects between a moving clock in 
a satellite and a stationary Time Receiver on the ground. This invention is a new and non- 
traditional method for performing a synchronization process that enables instantaneous 
synchronization of such Time Receivers (clocks) to a high level of precision. One use of 
this invention could be to modify an existing clock to perform the synchronization 
process defined in this patent application. 

As refered to above, the problem addressed in this disclosure is the instantaneous 
synchronization of two or more clocks located at different positions in a given spatial 
coordinate system, by means of either two-way or one-way radio communication 
between clocks. Generally, some or all of the clocks will be in motion relative to that 
coordinate system. In a typical application the coordinate system is attached to the earth, 
and some of the clocks would be on the surface of the earth, while others would be in 
satellites orbiting the earth. The satellite clocks would typically be atomic (Rubidium or 
Cesium), and the earthbound clocks might be either crystal oscillator or atomic. 

When some of the clocks are in moving satellites, in one aspect of the inventio, 
synchronization requires that relativistic effects be taken into account. 

Einstein's theory of Special Relativity is based on two postulates. The first is 
commonly called the relativity postulate, which states that for any two inertial frames in 
uniform relative motion, the laws of physics within each frame are identical. The second 
postulate is that the speed of light is constant, regardless of the inertial frame in which it 
is measured. 
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Einstein showed that if both postulates are assumed, that neither time nor distance 
is absolute. Specifically, if measurements of the time interval and/or distance between 
two events in space-time are made within each of two inertial frames in relative motion, 
the measurements will differ. 

The non-absolute nature of time posited by the Special Theory of Relativity 
poses complications in the process of time synchronization, and has led to errors resulting 
from misapplication of the Special Theory. 

The invention described herein removes the complications imposed by the Special 
Theory, by assuming only the first of Einstein's two postulates to produce a space-time 
theory in which time can be regarded as absolute. As a consequence, distance and 
velocity are defined somewhat differently, but the process of time synchronization is 
made extremely simple. 

SUMMARY OF THE INVENTION 

In one aspect the invention is a method for synchronizing one clock to another 
clock or a plurality of clocks to a given clock; in one particular case where one clock 
defined as a reference clock is on the earth and another clock is in a satellite which is in 
translatory motion relative to the reference clock. In that case, the method is to 
synchronize the satellite clock to the reference by making measurements of times of 
transmissions and arrival of a sequence of two-way transmissions and using those 
measurements deriving a correction term to be applied to the satellite clock. In one 
aspect of this case, the satellite clock is reset ascending to the correction and in another 
aspect of this case the correction term is used along with the satellite clock reading. 
Further, in the latter aspect the correction term may be applied in the satellite and a 
corrective time broadcast or in an alternative the satellite time and the correction term are 
both broadcast to user equipment which performs the correction. 
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Another primary aspect of the invention is referenced to as one-way 
communication, meaning communication from a satellite to two or more ground stations 
whose clocks are designed to be synchronized. For convenience, one of the ground 
stations is designated as the reference station or clock. Notably, where there is a plurality 
of ground stations, they can all be synchronized to the reference station simultaneously, 
or by designated groups or serially, depending on the communication system between the 
reference station and the other stations. In the one-way communication aspect first a 
corrective term representing difference in time reading between the reference station and 
the satellite is obtained, then a corrective term representing the difference in time reading 
between the next station, or stations and the satellite is obtained, then the difference 
between the corrective term is obtained, one that difference is applied to the station 
(clock) to be synchronized to the reference station. 

In another aspect, applied to a specific system such as GPS, the invention is 
implemented between a GPS ground control station and the satellite to determine a 
correction term for each satellite. Then the satellite can use the corrective term locally (at 
the satellite or broadcast it along with a time signal so that the correction calculation is 
done by the user). 

BRIEF DESCRIPTION OF DRAWINGS 

Figure 1 SpaceTime diagram for synchronization using two-way communication 
Figure 2 Correction (synchronization) of satellite clock using two-way 



Figure 3 



communication and subsequent synchronization of a clock C 
SpaceTime diagram for synchronization of clocks using one-way 



communication 



Figure 4 
Figure 5 
Figure 6 



Synchronization of clocks using one-way communication 
Flowchart for synchronization using one-way communications 
Apparatus for synchronization of clocks using one-way communication 
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DETAILED DESCRIPTION 

Following are descriptions of methods and apparatus for synchronizing clocks 
according to the invention. 

In describing the synchronization process, it is convenient to regard one of the 
clocks as a reference to which other clocks will be synchronized. This reference clock 
will be called clock A in the following discussion. In some applications, clock A is 
located on the surface of the earth, and the clocks to be synchronized are either on the 
surface of the earth or in satellites orbiting the earth. Also clock A could be located in an 
orbiting satellite. However, for convenience of the explanation, clock A is described as 
being located on the surface of the earth. 

If another clock is to be synchronized to clock A, and there are no other clocks 
available, then two-way line-of-sight radio transmissions may be employed to 
synchronize the second clock, but one-way transmissions are not enough to accomplish 
this unless additional information is available. It is also necessary that the clocks be 
visible to each other. On the other hand, it may be possible for the second clock to be 
synchronized with only one-way transmissions from a third clock which has both clock A 
and the second clock in view, even if the latter two clocks are not be visible to each other. 
This will generally require that the third clock broadcast certain information regarding its 
position and events (epochs) of signal transmission as a function of time. Such is the 
case if the third clock is located in a Global Positioning System (GPS) satellite or a 
GLONASS satellite. 

Synchronization Using Two-Way Radio Transmissions 

Understanding of the following explanation will be aided by referring to Figures 1 
and 2. Consider the case involving three clocks, denoted by A, B, and C. Clock A is 
considered to be the reference clock, and clocks B and/or C are to be synchronized to 
clock A. Clocks A and C are visible from clock B, but clocks A and C may not be visible 
to each other. Although clocks B and C might be moving relative to clock A, we assume 
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here that clock C is fixed relative to clock A, and that clock B has uniform motion 
relative to the other two clocks. 

Co-located with each clock is a radio station and a computer. It is assumed that 
the radio station in clock A is capable of two-way communication with the radio station 
in clock B, and that clock B can send a one-way signal to clock C. For mutational 
simplicity, the letters A, B, and C will also be used to identify the locations of the clocks, 
as well as the radio stations/computers associated with each clock. 

Clocks A and C are assumed to be at rest with precisely known locations within a 
3-dimensional inertial (i.e. non-accelerating) coordinate system, and clock B has uniform 
translatory motion (i.e., constant velocity) relative to the inertial system. All three clocks 
are identical in the sense that they would run at precisely the same rate if they were all 
stationary within the same inertial frame. However, the clocks may have been initialized 
differently, so their times may not agree. 

We first describe the two-way synchronization of clock B to clock A followed by 
a synchronization of clock C to clock B, from which a synchronization of clock A to 
clock C can be established. The synchronization process uses absolute time within the 
new space-time model described above. 

Measurement of Relative Motion 

Before synchronization can be achieved, measurements must be made which 
characterize the radial velocity of clock B relative to clock A. These measurements are 
made using clocks A and B, and do not involve clock C. An important parameter needed 
for the synchronization process is called the characteristic parameter of relative motion, 
and is denoted by This parameter is a measure of the relative radial velocity of clocks 
A and B, which is the rate at which the distance between A and B is changing with 
respect to time. To define % 9 assume that A sends two signals to B which are spaced At A 
seconds apart according to clock A. Due to the relative motion of clocks A and B, the 
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two signals will generally arrive at clock B with a different time spacing as measured by 
clock B. The parameter £ is simply the ratio of the latter time spacing to the former, i.e., 
the two signals arrive with time spacing ^At A according to clock B. Because the relative 
motion of A and B is uniform, £ does not depend on the times that the two signals were 
transmitted from A; in particular £ does not depend on At A . Note that if there is no 
relative motion between A and B, then the time spacing between the two signal arrival 
times is equal to the spacing between the transmission times; in this case, £=1 . 

On the other hand, suppose that the direction of signal transmissions is reversed, 
so that B sends two signals to A which are spaced At B seconds apart according to clock B. 
According to the relativity postulate, the symmetry of relative motion of A and B will 
result in the two signals arriving at clock A with time spacing ^A B according to clock A. 

The parameter % can be determined from a sequence of measurements made by 
sending radio signals back and forth between A and B. This process is diagrammed in 
Figure 1 . Suppose a radio signal S w is transmitted from A at time tf according to clock 
A, with the time of transmission tf as a message for B. The signal S w is received at B at 
time t? according to clock B. B instantly records tf and rf , and simultaneously 
transmits back to A a signal S x containing time t? as a message for A. The signal S x is 
received at A at time t* according to clock A. A instantly records and , and 
simultaneously transmits back to B a signal S y containing time as a message for B. 
The signal S y is received at B at time t\ according to clock B. B instantly records t£ and 
t\ , and simultaneously transmits back to A a signal S z containing t\ as a message for A. 
This signal arrives at A at time tf according to clock A. At this point in time, A has 
knowledge of tf , t£ , t* , ff , and t\ , which constitute all the basic measurements 
needed to establish the relative motion parameter %. It should be noted that B also knows 
all of these measurements, except for t* , and could have knowledge of t* by an 
additional message transmitted from A. 
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From the definition of £ given above, we see that 



l 2 l x l 2 t x 



(1) 



Since depends only on the ratio of time differences, it is independent of how 
the clocks were initialized, so that it can be determined before any attempt has been made 
to synchronize clock B to clock A. Taking the product of the two expressions for ^ above 
and taking the square root, we can obtain an expression for ^ in terms of only the 
reception times at clock A: 



l _3 h_ 

V2 h J 



(2) 



Definition of Synchronization 



The process of determining ^ normally can be done so rapidly that the relative 
radial velocity of clocks A and B can be considered constant during the interval that the 
measurements are being taken. If the relative radial velocity of clocks A and B 
thereafter remained constant and nonzero, at some point in time (possibly prior to the 
measurements described above) the distance between A and B would be zero, i.e., the 
clocks would coincide. Such an event is fictional, because the trajectories of the clocks 
will never permit them to be at the same position at any point in time. However, by 
assuming a constant relative radial velocity we can define the fictional moment of spatial 
coincidence which is useful in defining synchronization. We define clocks A and B to be 
synchronized if they would indicate exactly the same time at this fictional moment of 
coincidence (the situation when A and B have no relative motion will be dealt with later). 

In the above sequence of signal transmissions and receptions, the times ff and t* 
are measurements made by clock B before it has been synchronized. If B were 
synchronized to A, the time readings /f and t\ would become sf and s\ ("s" is used to 
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denote "synchronized"), but tf , , and ^ would remain the same. Clock B can be 
synchronized by resetting it so that it would read af at the moment it would otherwise 
read ff . This is accomplished by determining , which determines the correction 
sf -ff that needs to be applied, defined as e B . 

One determines s* by assuming the clocks were synchronized, so that each 
would indicate the same time f 0 at the fictional moment of spatial coincidence. Imagine 
that A sends a radio signal at that very moment. The signal is simultaneously received at 
B at time t 0 according to synchronized clock B. Let us also consider the signal sent by A 

at time ff and received at unsynchronized time ff , as previously discussed. This 
reception time is sf according to synchronized clock B. Using the definition of % 9 we 
have 



On the other hand, still assuming the clocks synchronized, imagine that at the 
moment of coincidence, B sends a radio signal at time f 0 according to clock B. The 

signal is simultaneously received at A at time f 0 according to clock A. Let us also 

consider the signal sent by B at unsynchronized time ff and received by A at time t£ , as 

previously discussed. The transmission time is sf according to synchronized clock B. 

Again, according to the definition of we have 



(3) 




(4) 




Taking the product of expressions (3) and (4), we obtain 
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(5) 



Then, 



1 



{?tt-4). 



(6) 



Substituting (6) back into (3) and performing algebraic simplification results in a 
solution for the value that jf that ff should have if clocks A and B were synchronized: 



Note that if there is no relative motion between A and B, then expression (7) will 
still give a correct value for s* by using £ = 1 , even though there may be no point in time 
where clocks A and B coincide. 

It can now be seen that as soon as tf , ff , and ^ and t\ have been measured, the 
parameter \ can be computed from (1), and then s* can be computed from (7). 
Synchronization of clock B can be accomplished in two equivalent ways: 

Method 1: Add the correction term s B = sf -ff to any unsynchronized reading of 
clock B. 

Method 2: Reset clock B by incrementing its indicated time by e B seconds. 

After synchronization is achieved, we now regard time as absolute, i.e., we freely 
use times observed by clocks A and B to calculate distances and velocities in the same 




(7) 
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way that would be done in a Newtonian system, that is, one in which relativistic time and 
distance dilation is ignored. 

At this point we assume clock B has been synchronized by method 2 above. After 
the clock is reset, its time readings will be denoted with s rather than t. For example, 
clock B's reading of t* at the instant of reception of the signal S w sent by A is replaced 
by s* . Without loss of generality, we may also assume that t 0 = 0 in expression (6), 
since the value of t 0 never appears in expression (7). 

Define the starred distance df B * from A to B at the instant sf of reception of 
the signal sent by A at time tf , as follows: 



where c is the velocity of light in the 3 -dimensional coordinate system as it travels 
from A to B. 

After substituting expression (7) for sf into (8), we obtain the alternate 
expression 




(8) 




(9) 



Now define the starred radial velocity vf 



,AB* 



between A and B as follows: 




(10) 
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An alternate expression for vf* can be obtained by substituting t Q = 0 in (5) to 

obtain 



4=<f> en) 



and then substituting into (10) to obtain 



vf=c^. (12) 

With starred radial velocity defined, the starred distance from A to B at an 
10 arbitrary time t B according to clock B is given by 

d A» =v A» s B 9 (13) 

where uniform translatory motion is assumed. 

15 

It should be noted that the above definitions of starred position and starred radial 
velocity are different from those that would normally be used within an inertial frame. 
An important consequence of these nontraditional definitions is that the speed of light in 
propagating from A to B (uplink speed) has the traditional value c, but the speed from B 
20 to A (downlink speed) is c- vf* . This can be shown as follows: Treating time as 
absolute, the downlink speed of light c is 



» _ distance between A and B at time sf _ df B * 
propagation time from B to A t 2 -s x 



25 Substitution of expression (8) for yields 



s: B -t A 

c *= c irA- (15) 
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Substitution of expression (7) for jf , followed by algebraic simplification, gives 

* c , 

o (16) 

By solving expression (12) for £ in terms of vf B * and substituting into (16), the final 
result is 

c=c-vf\ (17) 

This result is important in the synchronization of clock C, which will now be described. 
It will also be used in the method of synchronization using one-way signal transmissions, 
to be described later. 



Using clock B to synchronize clock C with clock A 

Clock C can now be synchronized to clock A by receiving a one-way downlink 
signal transmission from B with a message containing the transmission time s B 
according to synchronized clock B, and ephemeris data which permits determination of 
the velocity vector v B and the position of clock B at time s B in the coordinate frame 
containing clocks A and C. Since the position of clock C is precisely known, the radial 
velocity \ CB and starred distance d CB * from C to B can be determined at time s B . Let 

the arrival time of the signal as measured by unsynchronized clock C be denoted by f . 
The signal arrival time s c measured when clock C is properly synchronized is given by 

d CB * 

Clock C may now be synchronized to clock B (hence to clock A) by one of two methods: 
Method 1: Add the correction s c = s c -f to any unsynchronized time of clock 

C. 

Method 2: Reset clock C by incrementing its indicated time by s c seconds. 
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Synchronization Using One- Way Signal Transmissions 

5 Understanding of the following explanation will be aided by reference to Figures 

3, 4, 5 and 6. Synchronization of two clocks using only one-way signal transmission 
from a third clock is also possible using the concepts of this invention. Here the two 
clocks to be synchronized will be denoted by A and C, and the third clock by B. A radio 
receiver is co-located with clock A and a radio receiver is also co-located with clock C. 
1 0 Co-located with clock B is a transmitter. 

1^ For simplicity of explanation, clocks A and C are assumed to be at rest with 

|3 precisely known locations within a three-dimensional coordinate system fixed to the 

|t! earth, and clock B may be moving within this system. A typical coordinate system is the 

1 1 1 5 GPS Earth-Centered, Earth Fixed (ECEF) system with orthogonal x- 9 y- 9 and z- axes 

■IjJ meeting at an origin located at earth's center. The coordinates of clocks A, B, and C will 

k3 

* be denoted respectively by (x A 9 y A 9 z A ) 9 (x\y B ,z B ), and (x c ,j/ c ,z c ). 

is 
ill 
IU 

%j Clock B broadcasts a radio signal which can be received by clocks A and C. The 

Sj 20 signal contains a navigation message with satellite ephemeris data which permits A and C 
to know the precise location (x 5 , y B , z B ) and velocity vector v B of clock B within the 

ECEF coordinate system at every moment of time as read by clock B. The signal also 
contains periodically occurring identifiable points in the signal structure called epochs, 
and the navigation message contains information which permits clocks A and C to 
25 identify the precise transmission time of each epoch according to clock B. 

Synchronization of clock C with clock A can be accomplished in a procedure that, 
for convenience, can be defined as having three stages. 

Stage 1 could be defined as a process for determining a correction for clock A that 
30 would be needed at clock A as if clock A was to be synchronized to satellite clock B. 
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That correction is defined as S . Stage 1 consists of calculations which normally take 

place at clock A. Figure 6 is a block diagram of an apparatus for performing the 
calculations of stage 1 . 

Stage 2 could be defined as a process for determining a correction for clock C that 
would be needed at clock C as if clock C was to be synchronized with satellite clock B. 

That correction is defined as S . Stage 2 consists of calculations that are similar to those 

of stage 1 and which normally take place at clock B and stage 2 uses an apparatus 
identical to that used for stage 1 for the stage 2 calculations, so Figure 6 can also be used 
as a block diagram for the stage 2 calculations by changing every occurrence of the 
superscript A to superscript B. 

Stage 3 defines a process of determining a correction term to be applied to clock 

CA 

C to directly synchronize it with clock A. That correction term is S 

Stage 1: Clock A receives the signal broadcast from clock B and records the time 
ff of reception of a specific epoch according to clock A. Clock A also extracts from the 

navigation message and records the time ff of transmission of the epoch according to 
clock B. It should be noted that time as indicated by clock B is proper time, that is, it is 
uncorrected for relativistic effects, and does not need to be synchronized with clock A. 
Clock A uses the ephemeris data contained in the received signal to calculate the precise 

location (xf,>f ,zf ) and velocity vector vf = (vf ,vf y ,vf ) of B at the moment ff of 
epoch transmission. 

At clock A the distance d M from clock A to clock B at the moment ff of epoch 
transmission is calculated: 
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Also calculated is the radial component vf of relative velocity between A and B 
(the component which lies along the line of sight from A to B) according to the formula 

v f=vf.u A \ (20) 

5 

where u M is a unit vector pointing along the line of sight from A to B at the 
moment of epoch transmission, and the dot indicates the dot (or inner) product of the two 
vectors. A positive value of vf corresponds to a distance from A to B which is 
increasing with time. An alternative method of determining vf is to compute it by 
1 0 measuring the signal Doppler shift Af in Hertz at clock A and using the formula 

I Af > (21) 

111 " " ^ 



W 



III 



= -c- 

f 



$i where /is the carrier frequency transmitted by the satellite, X is the corresponding 

11 1 5 wavelength at the propagation speed c of light in free space. A positive value of 
J/1 Af corresponds to negative value of vf B . 



At clock A the propagation time of the signal epoch in traveling from B to A is 
calculated according to the formula 

20 C=tV- (22) 

v v r 



This formula is different from common practice in that the denominator is not the 
speed of light, but the speed of light reduced by the relative radial velocity of the satellite. 
The reason for this difference is that, in common practice, an attempt is made to correct 
25 the readings of clock B for relativistic effects of satellite motion and earth rotation 
relative to an inertial frame, prior to computing the propagation time t B p f op . However, 
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such corrections are unnecessary and lead to misapplication of the Lorentz 
transformations developed by Einstein. By using only proper times (i.e., uncorrected 
times of freely running clocks) and regarding time as absolute, no corrections of the 
satellite clock is necessary, but the formula for t BA op changes to that shown above. 

5 An epoch arrival time sf at clock A is now calculated, which is based on clock 

B's transmission time and the signal propagation time, in contrast to the direct 
measurement of arrival time t A by clock A: 



§ 10 

.1*1 



m 



mi 



III 



•*1 l \ ^ l prop 



B d AB • (23) 



1 c-vf 



The last step in the calculation at clock A is to compute the correction term 

e A =sf-tf (24) 

for clock A. At this point, clock A could be synchronized to clock B in the 
15 satellite if desired, either by adding s A to any unsynchronized reading of clock A, or by 



Hi resetting clock A by incrementing its reading by s A as soon as s A has been determined. 

Q However, if clock A is itself a time reference station, such a synchronization would not 



normally be performed. 



20 Stage 2: Stage 2 is identical in all respects to the calculations made by clock A, 

except that the calculations use data from the signal as received at clock C. Figure 6 
serves as a block diagram of an apparatus for performing these calculations if every 
occurrence of superscript A is replaced with superscript C. 

25 Clock C receives the signal broadcast from clock B and records the time t% of 

reception of a specific epoch according to clock C. This epoch may be the same epoch 
record by clock A, or it may be a different epoch. Also the epoch could be transmitted 
before or after the epoch used by clock A. Clock C also extracts from the navigation 
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10 



message and records the time t B of transmission of the epoch according to clock B. 
Clock C uses the ephemeris data contained in the received signal to calculate the precise 
location (x B ,y B ,z B ) and velocity vector v B = (v% x ,v% y ,v B z ) of B at the moment t B of 
epoch transmission. 

At clock C the distance d CB from clock C to clock B at the moment t B of epoch 
transmission is calculated: 



d CB - i4-* c h{yl-y c h{4-z c ) 2 , (25) 



Also calculated is the radial component v r of relative velocity between clock C 
and clock B according to the formula 



f v c r B =v**ii c \ (26) 



1 5 where u CB is a unit vector pointing along the line of sight from clock C to clock B 

at the moment of epoch transmission, and the dot indicates the dot (or inner) product of 
W the two vectors.. A positive value of v c r B corresponds to a distance from clock C to clock 

13 B which is increasing with time. An alternative method of determining vf* is to 

compute it by measuring the signal Doppler shift Af in Hertz at clock C and using the 
20 formula 

__M , (27) 

/ 

where/is the carrier frequency transmitted by the satellite, X is the corresponding 
25 wavelength at the propagation speed c of light in free space. A positive value of 
Af corresponds to negative value of v c r 
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25 



At clock C the propagation time of the signal epoch in traveling from clock B to 
clock C is calculated according to the formula 



d CB 

t pwp= CB ' ( 28 ) 

c-vjr 



An epoch arrival time ^ at clock C is now calculated, which is based on clock 
B's transmission time t\ and the signal propagation time, in contrast to the direct 
measurement of arrival time t% by clock C: 



C +B . +BC 



prop 



W D d CB . (29) 

iu 15 



The last step in clock C's calculations is to compute the correction term 

s c =s c 2 -t% (30) 

20 for clock C. 

Stage 3: At this point, clock C could be synchronized to clock B in the satellite if 
desired, either by adding e c to any unsynchronized reading of clock C, or resetting clock 
C by incrementing its reading by s c as soon as s c has been determined. 



In order to synchronize clock C to clock A, the correction term s is computed 
as follows: 
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6 c *=e c -s A . (31) 



The synchronization may be performed adding s CA to any unsynchronized 
reading of clock C, t c , or resetting clock C from t c by incrementing its reading by s CA 
after s CA has been determined. Normally the calculation of s CA is performed at clock C; 
in that case the value of s A (or the measurements at clock A required to calculate it) 
would be communicated from clock A to clock C by an independent communications 
link. 

Therefore the synchronized time of clock C to clock A is given by: 

S c = t c + S CA (32) 



In other words, the synchronization use of 8 may be implemented in either of 2 

ways: 

1 . Clock C may simply be reset according to S CA ; or 

2. Clock C is allowed to continue its prior readings which are modified by 
adding S CA to its reading and using the result as a synchronized time from 
clock C. 

For example, in the GPS, the satellite broadcasts the set reading of its clock plus a 
correction term and the user receiver applies the correction term to the set reading to 
achieve a corrected time. 



Description of Figure 6 



Figure 6 shows an apparatus for determining the correction term s A which is used 
in the calculation of the correction term s CA . An second identical apparatus is used to 
determine the other correction term s c needed for the calculation of s CA . Thus, Figure 6 
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also describes the second apparatus if all occurrences of superscript A are changed to B, 
and all occurrences of subscript 1 are changed to 2. 

An antenna (10) is connected to a satellite signal receiver (12) which receives 
signals from a satellite containing clock B. Clock A (14) is a continuously running clock 

which provides a continuous time reading t A which is available to the satellite signal 
receiver (12). The output t A of clock A (14) may also be fed to an optional time corrector 
(16). Clock A (14) could be an atomic clock or a clock governed by a crystal oscillator. 
The satellite signal receiver (12) provides three types of data, ff , t A , and satellite 
ephemeris data, which are stored in computer memory (18). Also stored in computer 
memory (18) is the known position [x A ,y A ,z A ) of the satellite signal receiver (12). 

The data tf , t A , satellite ephemeris, and {x A ,y A ,z A } in memory are available to 

a central processing unit (20), which has software modules (22), (24), (26), and (28) used 
in computing s A . The time t* is the transmission time of an identifiable epoch 
contained in the signal from the satellite containing clock B, and time t A is the value of 

the time reading t A from clock A (14) at the moment that the epoch is received. The 
satellite ephemeris data is information which permits the calculation of the position and 
velocity of the satellite at any given time. 

Module (22) is a satellite position and velocity calculation in which t* and the 
satellite ephemeris data are fetched from memory (18) to produce the satellite position 
( x i>yi> z \ ) velocity vector vf at time t* . The position and velocity vector are fed 
to module (24), which is a distance and radial velocity calculation. This calculation also 
requires the known position {x A ,y A ,z A ) of the satellite receiver (12) which is fetched 

from memory (18). The outputs of module (24) are the distance d AB of the satellite 
receiver (12) from the satellite containing clock B and the radial velocity v AB of the 
satellite containing clock B at time *f . These outputs are fed to module (26) which is a 
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calculation of the propagation time t prop of the signal epoch as it travels from the satellite 
to the satellite signal receiver (12). The propagation time t p A p is fed to module (28), 
which is a calculation of s A . The calculation of e A also requires ff and t A , which are 
fetched from memory (18). The output s A from module (28) is stored in memory (18). 
The optional time corrector (16) can fetch s A from memory (18) to correct clock A (14). 

Also, s A is available for the subsequent calculation of e CA which is used to synchronize 
clock A (14) to clock C. 

Although particular embodiments of the invention have been described and 
illustrated herein, it is recognized that modifications and variations may readily occur to 
those skilled in the art, and consequently it is intended that the claims be interpreted to 
cover such modifications and equivalents. 
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